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ABSTRACT 

We study the large N degeneracy in the structure of the four-point amplitudes of i- 
BPS operators of arbitrary weight k in perturbative M = 4 SYM theory. At one loop 
(order g 2 ) this degeneracy manifests itself in a smaller number of independent conformal 
invariant functions describing the amplitude, compared to AdSs supergravity results. To 
study this phenomenon at the two-loop level (order g 4 ) we consider a particular M = 2 
hypermultiplet projection of the general M = 4 amplitude. Using the formalism of N = 2 
harmonic superspace we then explicitly compute this four-point correlator at two loops and 
identify the corresponding conformal invariant functions. 

In the cases of ^-BPS operators of weight k = 3 and k = 4 the one-loop large N 
degeneracy is lifted by the two-loop corrections. However, for weight k > 4 the degeneracy 
is still there at the two-loop level. This behavior suggests that for a given weight k the 
degeneracy will be removed if perturbative corrections of sufficiently high order are taken 
into account. These results are in accord with the AdS/CFT duality conjecture. 
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1 Introduction 



In recent years the supersymmetric quantum conformal field theories in dimensions higher 
than two have undergone a renaissance due to the discovery of the AdS/CFT correspon- 
dence [1]. In the first place this concerns N = 4 super Yang- Mills (SYM) theory in four 
dimensions whose conjectured dual is type IIB (superstring) supergravity on an anti-de 
Sitter background. 

On the string theory side the present investigations are focused either on low-energy 
supergravity or on semiclassical string solutions carrying non-vanishing global charges of 
the superconformal group. In the latter case a direct comparison of the corresponding string 
energies to the spectrum of scaling dimensions of certain composite operators in N = 4 SYM 
is possible [2, 3, 4]. 

According to the duality conjecture, AdS supergravity is dual to the gauge theory taken 
at an infinite value of the 't Hooft coupling A and in the large iV limit, where iV is the rank 
of the gauge group (S)U(JV). In this limit the gauge theory reveals the following peculiar 
behavior: 

• the sector of protected operators is enlarged; 

• the spectrum of anomalous dimensions degenerates; 

• the symmetry is enhanced. 

The first property is a simple consequence of the rather general fact that in the large N 
limit the anomalous dimension of a multi-trace operator made of single-trace constituents 
is equal to the sum of their individual anomalous dimensions. In this sense, in the large N 
limit the protected operators form a ring 1 . 

The second property is more subtle and it is reminiscent of the free theory pattern, where 
usually several operators with the same canonical dimension coexist. It is well known that 
protected (for instance, ^-BPS) operators with the same quantum numbers can be realized 
in the gauge theory both as single- and multi-trace operators. Thus, multiplying protected 
operators with the same quantum numbers by an unprotected operator, we obtain different 
operators with the same dimension in the planar limit, hence the degenerate spectrum. 
Such a mechanism is essentially due to the conformal supersymmetry responsible for the 
protection and, again, due to the large N factorization property of the Feynman graphs. 
However, this is not the whole story. As was recently shown in Ref. [5], there exist, for 
instance, single-trace operators (which cannot be obtained as the product of a protected and 
an unprotected operators) whose anomalous dimensions coincide up to two loops. This is 
evidence for an extra symmetry in the planar limit which is responsible for this degeneracy 
of the perturbative spectrum. The symmetry enhancement [5] is an interesting property of 
the gauge theory related, for instance, to the integrability of the spin chain Hamiltonian 
[6]. The latter can be identified with the one- loop dilatation operator in N = 4 SYM. 

Recently, another interesting peculiarity of the large ./V limit of N = 4 SYM was found 
[7]. We recall that in regard to the supergravity approximation of the gauge theory, the 

1 When computing the two-point function of a multi-trace operator, one can readily see that the leading 
1/N contribution is due to the Feynman graphs which contribute separately to the two-point functions of 
the individual single-trace constituents. 
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most interesting object to study are the four-point correlation functions of ^-BPS operators 
(for reviews see, e.g. [8, 9]). These operators are dual to the Kaluza-Klein states of the 
compactified type IIB supergravity. They form isolated series of superconformal representa- 
tions and their scaling dimension as well as their two- and three-point correlation functions 
are protected from receiving quantum corrections. In contrast, their four-point correlators 
undergo perturbative renormalization. 

The four-point amplitude of ^-BPS operators with arbitrary weights (dimensions) k p , 
p = 1,2,3,4 is specified by a number of conformal invariant functions depending on the two 
conformal cross-ratios. There is a general procedure based on the field-theoretic insertion 
formula which allows one to determine the maximal number of a priori independent such 
functions appearing in the quantum part of the amplitude [10, 11]. 

The explicit amplitudes corresponding to the cases where all k p = k are equal to 2, 
3 or 4 have been computed in the supergravity approximation by using the effective five- 
dimensional description of type IIB, and also in perturbation theory at one loop (order g 2 ) 
and in the planar limit. It was then observed in Ref. [7] that the perturbative amplitudes 
exhibit a certain degeneracy in comparison to their supergravity counterparts. Starting 
from k = 4, the number of conformal invariant functions describing the amplitude at one 
loop is smaller than in the supergravity regime. We can formulate this fact as another 
interesting property of the large N limit: 

• the perturbative four-point amplitudes of ^-BPS operators degenerate. 

Thus, not only the two-point functions (the spectrum) but also the higher-point perturbative 
amplitudes exhibit simplified features in the large N limit. 

The aim of the present paper is to the study this degeneracy in more detail. Recall 
that for the first time the degenerate situation occurs for k = 4. In the supergravity regime 
one finds two different conformal invariant functions, precisely the number predicted on 
general grounds. However, at one loop the two functions describing the large N amplitude 
coincide. It should be noted that the perturbative computation has been performed for ^- 
BPS operators realized as single-trace operators. Generically, the ^-BPS operators can be 
realized as mixtures of single- and multi-trace operators with arbitrary mixing coefficients. 
Therefore, two natural questions can be asked about this large N degeneracy: 

1. Does the operator mixing influence the degeneracy? 

2. Is the degeneracy lifted when the higher-loop corrections are taken into account? 

We study the first question by using the k = 4 amplitude as the simplest example. The 
^-BPS operator of dimension 4 is a mixture of one single- and one double-trace operator, 
the latter is made of two single-trace constituents of dimension 2. Firstly, we show that 
the free amplitude involving mixed operators coincides with that for single traces in the 
large N limit, provided that the mixing parameter scales faster than l/yN. As discussed 
in Ref. [7], the supergravity-induced amplitude admits a unique splitting into a "constant" 
and an "interacting" parts. The former precisely matches the free amplitude computed for 
single-trace operators. In the AdS/CFT language the interacting part must comprise all 
loop corrections to the free amplitude. Clearly, in order not to affect this agreement with 
the constant part of the supergravity amplitude, in field theory we should allow only mixing 
which is more heavily suppressed than 1/y/N. Secondly, we argue that with such a mixing 
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the contribution of the double-trace operators at one loop is suppressed in the large N limit 
as well, so the degeneracy is not affected by the operator mixing. 

We then investigate the issue of the higher-loop corrections by extending the perturba- 
tive treatment to two loops (order g 4 ). To illustrate the more general picture, we consider 
the amplitude involving ^-BPS operators of equal but arbitrary weight k. By performing 
an explicit diagrammatic computation we find that at two loops the degeneracy occurring 
in the case k = 4 is lifted: Two different conformal invariant functions emerge. However, for 
amplitudes involving operators of weight k > 4 the degeneracy persists even at the two-loop 
level. These observations indicate that for a fixed weight k the degeneracy is removed when 
loop corrections of sufficiently high order (depending on k) are taken into account. Since 
the supergravity-induced amplitudes comprise all perturbative corrections, it is natural to 
expect them to realize the maximal possible number of different conformal invariants. 

In fact, the perturbative degeneracy phenomenon has a rather simple origin. At a given 
order of perturbation theory we have only a finite number of topologically different planar 
interacting graphs. Whether all of them are indeed present depends on the weight k of the 
four composite operators for which the amplitude is computed. It is clear that increasing k 
makes all possible interacting topologies appear. Eventually, we reach the point of saturation 
beyond which increasing the weight amounts to decorating the interacting graphs with free 
propagators. In the large N limit this does not generate new conformal invariants. On the 
contrary, fixing k and increasing the order of perturbation theory "awakes" new interacting 
topologies and hence generates new conformal invariants, thus lifting the degeneracy. 

One could also ask if there is any relationship between the degeneracy of the four-point 
amplitudes and that of the spectrum observed in Ref. [5] . Any conformal invariant function 
can be expanded over an infinite basis of conformal partial wave amplitudes, each of which 
represents the contributions of an individual primary operator. The degeneracy of the four- 
point functions is therefore related to some kind of degeneracy of the spectrum and of the 
OPE coefficients. It would be interesting to find out if the degeneracy of the spectrum is also 
removed when higher loop corrections are taken into account. Indeed, in the lattice picture 
[6] the L-loop dilatation matrix induces the mixing of the L+l neighboring elementary fields 
constituting a composite operator, i.e. with L growing the effective interaction is spreading 
over the lattice and can subsequently increase its size which is equal to the dimension of 
the operator. When this happens, the degeneracy might be lifted in a way similar to that 
for the four-point amplitudes. 

The organization of the paper is as follows. In Section 2 we describe the structure of the 
Af = 4 four-point amplitude of ^-BPS operators determined by the field-theoretic insertion 
procedure and discuss its "pure" N = 2 hypermultiplet projection. It involves a subset of 
k — 2 conformal invariant functions (or one such function for k = 2) from the original M = 4 
amplitude. For k > 4 this subset is not complete but it is sufficiently large to illustrate the 
degeneracy problem. The advantage of the hypermultiplet projection is that we can use the 
quantum formalism of M = 2 off-shell harmonic superspace [12] which is particularly efficient 
at two loops 2 . For the sake of clarity, at the end of Section 2 we present and discuss our main 
results, whose derivation is explained in detail in Sections 4 and 5. In Section 3 we discuss 
the influence of operator mixing on the degeneracy problem. In Section 4 we summarize the 
N = 2 insertion procedure and the Feynman diagram tools necessary to perform the two- 
correlation functions of i-BPS operators at order g 4 have also been computed in the Af = 1 superspace 
formalism [13, 14]. 



3 



loop computation of the hypermultiplet projection. In Section 5 a diagrammatic calculation 
is presented and the corresponding conformal invariants are identified. Finally, in Section 
6 we use the operator product expansion and the knowledge of the one-loop anomalous 
dimensions of certain operators to get an insights into the structure of the k = 3 and k = 4 
amplitudes, independently of the diagrammatics. In particular, for k = 3 we are able to 
completely reconstruct the two-loop amplitude from these OPE considerations. Appendix 
A clarifies the technical tool of harmonic analyticity which we systematically use in order 
to drastically reduce the number of relevant Feynman graphs. Appendix B contains the 
details of the graph calculation for k = 3 and 4 whose generalization for arbitrary k is given 
in Section 5. 



2 Generalities and main result 

The lowest component of a |-BPS multiplet in Af = 4 SYM is a real scalar field of dimension 
k transforming in the irrep [0, k, 0] of the R symmetry group SO (6) ~ SU(4:). In terms of 
the elementary fields it can be realized, e.g., as a single-trace operator 

Tr(^ ai . . . a »>) . (2.1) 

Here <j) a , a = 1, . . . , 6 are the N = 4 SYM scalars and {, } denotes traceless symmetrization. 
A convenient way to handle the SO(6) indices is to project the operator (2.1) onto the 
highest weight state of the irrep [0, k, 0]. This can be done with the help of a complex null 
vector u a (u a u a = 0) carrying U(l) charge, which can be viewed as a harmonic variable 
parametrizing the coset space SO(6)/SO(2) x SO(4) (see, e.g., Ref. [10] for details): 

O ik) =u ai ... u afc Tr(0 ai . . . <f) a ") . (2.2) 

We start by summarizing the general properties of the four-point amplitude of i-BPS 
operators in the N = 4 SYM theory. It is sufficient to restrict our discussion to the case 
where all the operators involved have equal weights k. On the basis of conformal and R 
symmetry alone the general form of the four-point amplitude can be parametrized as follows: 

(k\k\k\k)= y, a ™nM xm y nzl i ( 2 - 3 ) 

m+n+l=k 

where m,n,l are non-negative integers. This expression is a polynomial in the three ele- 
mentary propagator (Wick) contractions of the four points 

(12) (34) (13)(24) (14)(23) 

A. — g T~ ' * — 2 T~ ' ^ ~~ 2 T~ ■ ( ZA 

x 12 x 34 x 13 x 24 x 14 x 23 

Here (12) = (21) = u\u\, etc. denote the SO(6) invariant contractions of the harmonic 
variables at two different points. Every monomial in eq. (2.3) corresponds to a certain 
propagator structure built out of the propagators of the elementary fields (see Fig. 1). We 
refer to eq. (2.3) as to the representation of the amplitude in the propagator basis. The 
coefficients o!^ n i{s,t) are arbitrary functions of the two independent conformal invariant 
cross-ratios 

9 9 9 9 

_ x 12 ^34 f _ x 14 x 23 

b (y 9? ^ 9 9* 

x 13 x 24 x 13 x 24 



4 



Super confer mal symmetry puts additional kinematical restrictions on the coefficient 
functions in (2.3) which take the form of differential equations (superconformal Ward identi- 
ties) [15, 16]. Further, dynamical constraints on the amplitude are obtained through the pro- 
cedure of inserting the N = 2 SYM action. According to the "partial non-renormalization 
theorem" of Ref. [17], the quantum M = 2 amplitude takes a factorized form (see Section 
4 for details). Recently this insertion procedure has been generalized to the M = 4 case 
[10, 7]. The factorized form of the quantum corrections to the M = 4 amplitude is 

(*ifcifcifc>* Mat - = r n=a y, ^lM xm y nzl ^ ( 2 - 5 ) 

m+n+l=k— 2 

where R^f=4 is a universal polynomial prefactor carrying weight 2 at each point. Explicitly, 

Rm=a = s x 2 + y 2 + t z 2 + ( s -t-i)yz 

+ {l-s-t) XZ + {t-s-l) Xy. (2.6) 

All the dynamical information is thus encoded in the conformally invariant coefficient func- 
tions T^ nl {s,t) which are the subject of our subsequent investigation. 3 Substituting the 
universal prefactor into eq. (2.5) and comparing it with eq. (2.3), we read off the following 
expression for the quantum corrections to the coefficients a!^ni : 

+ (s-t-1) F {k) , , , + (l-s-t) F [k) , , , + (t - s - 1) F {k) . ... 

1 V / m,n— 1,1— 1 v ' m— l,n,l— 1 V / m— l,n— l,i 

(k) 

If any of the indices of in this equation becomes negative, then the corresponding term 
is absent. 

Since all the operators involved are identical, the four-point amplitude is invariant under 

the symmetric group S4 which permutes the points 1, . . . , 4. Only the subgroup 53 acts non- 

(k) 

trivially on the cross-ratios s and t and, as a result, on the coefficient functions F~,rd m ^ ^ s 
generated by the two independent crossing transformations, e.g., 1 2 and 1 <-> 3. If all 

(k) 

the indices m,n,l of F^nl are different, one can use the action of S3 to order them, e.g., 

(k) 

m > n > I. The function Fi^i, with m > n > I does not satisfy any crossing symmetry 
relation. It can be taken as a representative of the corresponding crossing equivalence class 

(k) 

consisting of six elements. If any two indices of J~~Li coincide, then the function satisfies 
an additional crossing symmetry relation (and the corresponding crossing equivalence class 
consists of three elements). Finally, in the case where all three indices coincide, the function 
transforms into itself under the whole group S3 (and the corresponding crossing equivalence 
class consists of a single element). 

3 To make the reader familiar with the notation we stress that the upper index in ^^„i(s,t) denotes the 
weight of the i-BPS operators. For the sake of clarity we sometimes separate the lower indices by commas. 
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Figure 1. The propagator structure corresponding to the monomial X m y n Z l . 

For sufficiently small values of k the number of conformal invariants J- representing the 
independent crossing equivalence classes is easy to work out (the general formula is given 
in Ref. [11]). For instance, both for k = 2, 3 the quantum part of the amplitude is given by 
a single function T ', while already for k = 4 two functions are needed. 

The functions T can be explicitly computed both in perturbation theory and in the 
strong coupling regime by using the AdS/CFT duality. In Ref. [7] it was observed that 
in the large N limit the perturbative four-point amplitudes of the ^-BPS operators exhibit 
a certain degenerate behavior in comparison to their supergravity partners: The number 
of independent conformal invariants describing the one-loop (order A ~ g 2 ) perturbative 
amplitude is smaller than that in the supergravity regime. In what follows we study this 
degeneracy phenomenon at the two-loop (order A 2 ~ g 4 ) level by explicitly working out the 
corresponding four-point amplitude. 

Computing four-point correlation functions at two loops is of course much more involved 
than at one loop, due to the rapidly increasing number of Feynman graphs. A very effi- 
cient technique is provided by the M = 2 harmonic superspace formalism. This approach, 
combined with the N = 2 insertion procedure, allows us to drastically reduce the number 
of relevant graph topologies to be computed. 

In order to apply the M = 2 harmonic superspace formalism we need to decompose 
the ^-BPS operators into their Af = 2 hypermultiplet (HM) and SYM constituents. This 
is done by first decomposing the N = 4 field-strength superfield W = u a (j) a {x) + ■ ■ ■ into 
the Grassmann analytic HM superfield q + = uf4> l (x) + • • • (and its harmonic conjugate 4 
q + = uf(p l (x) + • • • ) and the chiral field strength W = w(x) + • • • (and its antichiral 
conjugate W = w(x) + •••). Here the six real scalars 4> a (x) forming an SO(6) vector have 
been split into an SU(2) doublet 4> % (x), i = 1,2 and a complex singlet w(x). The SO(6) 
vector harmonics u a have been replaced by SU(2) fundamental harmonics (u~~ = (u + )* , 
ufe^uj = 1), parametrizing the sphere S 2 ~ SU(2)/U(1). 

In the sequel we will be interested in the "pure" HM projection in which we put only 
(ps at points 1 and 4 and only q's at points 2 and 3. As explained in Ref. [7], taking 
this projection means considering only the terms with I = in eq. (2.3) (in Fig. 1 they 
correspond to the graphs without diagonals). Further, the propagator structures (2.4) are 

4 Here q means the usual complex conjugation for the field together with an antipodal reflection on the 
sphere S 2 ~ SU(2)/U(1) for the harmonic variable. 
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replaced by 



X ^ X= BW*, y~Y=M&, (2.8) 

x 12 x 34 x 13 x 24 

where, e.g., 

[12] = -[21] ee {«i)t&(ua)t (2-9) 

is the SU(2) invariant contraction of the two sets of harmonic variables at points 1 and 2. 
Under this projection eq. (2.5) reduces to 

k-2 

(Q (*)) = (§V|g*|§*) = R M=2 ^\-m-2,M X m Y k ~ m ~ 2 , (2.10) 

m=0 

where 

R^=2 = sX 2 + {t- s-l) XY + Y 2 . (2.11) 

The factorized expression (2.10) with the universal prefactor (2.11) is precisely what one 
finds by the direct application of the field-theoretic insertion procedure to the four-point 
amplitude in the M = 2 theory (see Section 4) . 

According to the previous discussion of the crossing symmetry, we can take 3~mk-m-2 o 
with 



m_ = l/2(fc - 2) < m < k - 2 = m + for k even (2.12) 

and 



l/2(fe — 1) <m<k-2 = m + for k odd (2.13) 

as representatives of the different crossing equivalence classes of the coefficient functions of 
the four-point amplitude. 

It should be pointed out that for k > 4 the whole M = 4 amplitude cannot be restored 
from its pure M = 2 projection 5 and consequently, the functions J r ^\_ m _ 2 above do not 
form a complete set of a priori independent conformal invariants needed to parametrize the 
amplitude. However, the subset of conformal invariants emerging in the pure HM projection 
is sufficiently large to illustrate the large N degeneracy problem at the two-loop level. 

Among the chosen set of independent coefficient functions only one (for k odd) or two (for 
k even) admit a pair of coincident indices. The first such function (s,t) = J~k-2,a,o( s )t) 
is obtained by setting m = m+, so it exists for any value of fe. 6 It has the additional 
symmetry property 

F { +\s,t) = 1/t F ( +\s/t, l/t) . (2.14) 

5 To this end we would also need some of the "mixed" projections of the type q k ~ p q p . Obtaining such 
projections involves some amount of linear algebra; in addition, the corresponding two-loop graphs are more 
complicated. Therefore here we prefer to restrict ourselves to the "pure" sector. 

6 One could instead take J~^l_ 2 o w hich obeys ^' t '_2o( s i') = o(^> s )' ^he relation between the 

two functions is ^1_ 2 o o( s >^) = V s fc-2 o(V s ! V s )- 
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(k) 

The second function J-^^ m _ is obtained by setting m = m-, so it exists only for k even. 
In this case we have 

^!,™_, M) = V* ^2,m- 1 o(V«,*A) • (2-15) 
It is more convenient to introduce the functions 

F£\s, t) = ^L m „ 2i0 (t, s) , m + m + (2.16) 

since for m = m_ the function iF_ (s, i) = m _ (t, s) has the same crossing symmetry 
as JF^ : 

^i fc) (s,t) = 1/t ^\s/t,l/t). (2.17) 

The functions with m_ < m < m + do not a priori obey any crossing symmetry 

relations. Finally, we note that the case k = 2 is exceptional as all the three indices of J 7 ^ 
vanish. This function transforms covariantly under the whole group S3. 

In summary, the pure HM projection is parametrized by the chain of conformally invari- 
ant functions whose "boundaries" T± are subject to the crossing symmetry relation 

(2.17). The function J-^ exists only for k even. 

Concluding this section, we give the main result of our two-loop calculation, the details 
of which are presented in Sections 4 and 5. In the large N limit we find 

4*>0M) = ^{\{t + i)[* w (M)] 2 + <*> (2) M) + \^ 2 \s/tMt)} , (2.I8) 

where the range of m is given by Eqs. (2.12) and (2.13), in particular, for k odd the function 
J-^ is absent. In eqs. (2.18) A = g 2 N /4tt 2 is the 't Hooft coupling and the functions 
and $( 2 ) are the so-called one- and two-loop scalar box conformal integrals 

h =-^r* (1) (*'*>' ^ 

J x 15 x 25 x 35 x 45 x 13 x 24 

f xf 3 d 4 x 5 d 4 X6 {i^ 2 ) 2 (2) 

/ ry 2 2 2 ,-y> 2 ™2 ry, 2 ™ 2 2 ,yi 2 (?) \ / 

J x 15 x 25 x 35 x 56 x 16 x 36 x 46 x 13 x 24 

We see that in the large N limit all but one conformal invariants describing the pure HM 

(k) 

projection are equal. It is only the "boundary" function !F\ which is distinctly different 
from the others. 

Having found the pure HM part of the four-point amplitude of the ^-BPS operators, 
we can now understand what happens to the large N degeneracy at the two-loop level. 
According to Ref . [7] , at one loop and in the large N limit all conformal invariants emerging 
in the pure HM projection coincide: 

J%\s,t) =J%H*,t) = (2-21) 
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(3) 

This degeneracy manifests itself already in the case k = 3: The single function, !F\ , 
describing the one-loop amplitude obeys the extra crossing symmetry relation F\_ (s,t) = 

(3) 

•7+ {t, s) which is not required on general grounds. Just as in the one-loop case, the two-loop 

(3) 

function J-\' from eq. (2.18) obeys eq. (2.17) imposed by the symmetry of the amplitude, 
but the one-loop "bonus" crossing symmetry relation does not hold anymore. We therefore 
conclude that this symmetry of the k = 3 one-loop amplitude is destroyed by the two-loop 
quantum corrections. 

Another type of degeneracy appears in the case k = 4. On the general grounds of 
superconformal and crossing symmetry two a priori independent conformal invariants are 
needed in this case, while at one loop both of them appear to coincide in the large N limit 
[7]. According to eqs. (2.18), at the two-loop level we indeed find two conformally invariant 
functions J-±^ which are distinctly different. Again we observe that the one-loop degeneracy 
is lifted when the two-loop corrections are switched on. 

However, we also see that both types of large N degeneracy mentioned above still hold 
at the two-loop level for operators of higher weight, k > 4. Indeed, for k > 4 the functions 
J~m ^ coincide (they are equal to !F_ for k odd) and have the same bonus crossing symmetry 
(2.17) as F^. 

Finally, we mention that in the special case k = 2 we find 

•^M) = §2^ + t + l )\^ l) (^ t )] 2 ( 2 -22) 

+ (t/s, l/s) + $( 2 ) (S, t) + IfcP) (s/t, l/t) 

s t 
This result has been previously obtained in Refs. [18, 14]. 



3 Mixing and large N degeneracy 

In field theory |-BPS operators can be realized as mixtures of single- and multi-trace 
operators. Here we investigate whether the operator mixing can affect the degeneracy 
problem of four-point correlation functions. As an explicit example we consider the case of 
weight k = 4 operators. 

As discussed in Ref. [7], for ^-BPS operators of weight k < 4 the complete M = 4 four- 
point amplitude can be reconstructed from the "pure" N = 2 hypermultiplet projection, 
eq. (2.10). In particular, in the case k = 4 we can write down 

(Q (4) ) = £ aS 4 - mi oM)X m y 4 -"\ (3.1) 

m=0 

This projection involves five of the fifteen coefficient functions of the original M = 4 am- 
plitude, which belong to the three crossing equivalence classes (0040,0400), ( a i30' a 3io) an< ^ 

(4) 
a 220- 

In order to study possible mixing effects on we compute the coefficients at tree 
level. There, generically, the ^-BPS operator is of the form 

C~S + aD, (3.2) 
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where we have defined 



S = Tr(q 4 



D = Tr(q 2 )Tr(q 2 



(3.3) 



q being the N = 2 Grassmann analytic HM superfield and a an ^-dependent mixing 
parameter. 

We fix the normalization coefficient for O by requiring its two-point function to be nor- 
malized canonically. Computing separately the two-point functions (SS), (SD) and (DD), 
in the case of SU(N) color group we find (since we are interested in the color combinatorics 
we do not exhibit the space-time dependence of the correlators) 



(SS) = 4(N 2 - 1) Vi , 
(SD) = 8(N 2 -1)V 2 , 
(DD) = 8(N 2 -1)V3: 



where, for convenience, we have introduced 
Pi 

Therefore, the canonically normalized operator reads 



^(iV 4 -6iV 2 + 18), v 2 = ^(2N 2 -3), V 3 



(N 2 + 1) . 



O 



A(N 2 - 1) [Vx + 4aV 2 + 2a 2 P 3 ) ] 5 (S + aD) 



(3.4) 



(3.5) 



(3.6) 



The tree level contributions to the coefficients are given by diagrams with propagator 
structures as in Fig. 1, with / = 0, n = 4 — m, m = 0, • • • ,4. The coefficient functions 
a[,4Q and a^QQ correspond to disconnected diagrams built as product of 2-point function 
diagrams. Therefore, given the above normalization, they are trivially equal to 1. In order 
to determine the other coefficients we need to compute the actual four-point amplitudes 
for the mixed operator O. This can be achieved by evaluating separately the four-point 
functions (SSSS), (DDDD) and (SSSD), (SSDD), (SDDD) as well as the amplitudes 
obtained from the latter by permuting 5 and D. 

Since it is enough to compute one coefficient function within each crossing equivalence 
class, we concentrate, for instance, on the coefficient function a$ . It is convenient to first 
evaluate the building blocks of Fig. 2. 




These building blocks can be obtained from the tree diagrams of the 2-point amplitudes 
(3.4) by cutting a line. This operation simply amounts to inserting a 5 fll f )1 for the line 
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cut without affecting the color structure. Therefore, the three building blocks are still 
proportional to the polynomials given in eq. (3.5). Explicitly, 



SS 
SS 
SS 



2 4 Pi S aibl 
2 5 V 2 5 aibl 
2 5 P 3 S aibl 



(3.7) 



These building blocks are then combined to realize all possible configurations of single- and 
double-trace operators. The result is 



z (4) 

l 3l0 



J 130 



2 8 (iV 2 - 1) (Vl + 16a 2 P 2 2 + 4 « 4 ^I + 8aPiP 2 + lQa 3 V 2 V 3 + 4a 2 Pi7> 3 ) 



2 4(at2 _ !)2 + 4a7 7 2 + 2a 2 V 3 y 
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2 8 (iV 2 - 1) (Pi + 4aP 2 + 2a 2 V 3 ) 2 _ 
2 4 (iV 2 _ !)2 ( Pl + 4a p 2 + 2a 2 P 3 ) 2 ~ ^ 2 " 1 



(3.8) 



Being independent of the mixing parameter a, it gives the finite N result for the tree 
level part of the coefficient function a^ , independently of the presence of the double-trace 
operator. In the large N limit this expression agrees with the result of Ref. [7]. We note 
that the independence of a is due to the particular structure (3.7) of the building blocks. 



For k > 4 we observe the same phenomenon in the coefficients a 



(fc) 

(fe-l)lO' 



Next we concentrate on a 22 o- Again, the calculation can be carried out by first computing 
the building blocks depicted in Fig. 3. 






Fi gure 3. Building blocks for the coefficient o 22 o 



Introducing the concise notation Tr (T a T b ...T c T c i) = (ab...cd) for the trace of the color 
generators, the algebraic expressions for the building blocks of Fig. 3 can be written as 



SS 



SD 



(N 2 - 9) 



+ 

2 3 
+ 



N 

,27V 2 +9 
! iV 2 



[(ai<2 2 &i& 2 ) + (a 2 ai&i& 2 ) + (aia 2 b 2 bi) + (a 2 ai6 2 6i)] 
(aia 2 )(bib 2 ) + 4(ai6i)(a 2 6 2 ) + 4(ai6 2 )(a 2 &i) ; 



(aia 2 bib 2 ) + (a 2 ai&i6 2 ) + {a\a 2 b 2 bi) + (a 2 ai& 2 6i) + (aibia 2 b 2 ) + (a\b 2 a 2 bi) 
{2N 2 - 3) 



N 



{a\a 2 ){bxb 2 ) 



DD -» 2 3 [(N 2 + 3){a 1 a 2 )(b l b 2 ) + 2(a 1 b l )(a 2 b 2 ) + 2{a 1 b 2 )(a 2 b 1 )] . 
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All possible configurations of the four-point amplitudes of S and D are obtained by com- 
bining the building blocks of Fig. 3 in all possible ways. The last contraction between two 
blocks produces an extra factor of 4. Pulling out the common factor 2 s (N 2 — 1), we find 

(SSSS) -+ (^-y^ + 135-^ + ^) ss7>4, 
(SSSD) -» -^(lOiV 6 -54iV 4 + 216iV 2 - 405) =V 5 , 

(SSDD) -» -X; UN 6 + 6N 4 -63N 2 + 135) =V 6 , 
N z 

(SDDD) -> ^(2iV 6 + 57V 4 + 18iV 2 - 45) = P 7 , 

{DDDD} -► (A^ 6 + 5iV 4 + 19iV 2 + 15) = P 8 (3.9) 
The final answer for the coefficient a^o is then 



(4 ) 16 (P 4 + 4aP 5 + 6a 2 P 6 + 4a 3 P 7 + a 4 P 8 ) 
1220 ~ (N 2 -l)(Vi+4aV 2 + 2a 2 V 3 ) 2 



(3.10) 



We notice that by setting a = (i.e., by neglecting the double-trace contributions) we 
obtain the finite N result for the free amplitude involving only the single-trace operator: 

(4) single 16(iV 4 - f N 2 + 135 - g + f#) .„ 1 

° 220 " (7V2_i)(iV2- 6 + ^) 2 1 j 

In the large A/" limit this expression agrees with that in Ref. [7]. 

Next we analyze the large N limit of the general expression (3.10). Taking the large iV 
limit of Vi we obtain 

(4) 16 (iV 2 + 40(iVa) + 24(iVa) 2 + 8(Na) 3 + (iVa) 4 ) 
0220 " (N 2 + 8(Na) + 2(Na) 2 ) 2 ' (3 ' 12) 

If in the large limit the mixing parameter is chosen to be of order one, then the double- 
trace part of O provides a dominant contribution and one finds a^o ~ 1- 

We recall that according to our general treatment in Section 2 the coefficient function 
°22o( s ' admits a unique splitting on the non-trivial ("quantum") function of the conformal 

l 220 



cross-ratios constructed out of T± and an integration ("free field") constant al 



4J (s, t) = a { 2 % + F + (s, t) + F + (l/s, t/s) + (t - s - l)F-(t, s) . (3.13) 

The constant part of a^ (s,t) was computed in the supergravity regime [7] with the result 
a 22() SUgra = W> 1,e ' ^ scales as 1/N 2 . Therefore, the AdS/CFT correspondence does not 
allow the mixing parameter a to be of order one. If we pick a = k/^/N + . . ., then it will 

,( 4 ) 



affect the coefficient a 220 at order 1/N : 

,fi r(l + ^ 4 ). (3.14) 



/ 4 ) - iip ^ ^ 

l 220 - ^2 1 
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However, if we want a^o to match the large TV" value computed from supergravity exactly, 
we ought to require k = 0, i.e. the mixing parameter has to scale faster than 1/y/N. 
Moreover, we note that if a is expanded in the odd powers of 1/N, then the free amplitude 
has a 't Hooft type expansion in powers of 1/N 2 . Therefore, the simplest assumption for 
the mixing coefficient is a ~ 1/N + 




Figure 4. The leading large N contributions to the four-point amplitudes {SSSS) and 

(DDDD) at order A. 

Let us now briefly discuss the influence of the mixing on the large N degeneracy occurring 

2 

at the one- loop level (order A = |^). As shown in Ref. [7], harmonic analyticity, i.e. 
the superconformal kinematics, allows one to reduce the one-loop calculation of the four- 
point amplitude just to a single Feynman graph given in Fig. 4. This statement holds 
irrespectively of the intrinsic trace structure of the BPS operators and therefore it can be 
applied to the operator O. 

The various contributions to the coefficient a^ from the amplitude at order A are 
then obtained by specifying the operator (S or D) at each corner of the graph in Fig. 4. 
Using a double line notation for the propagators in there, the leading large N behavior is 
easily identified: Both the purely single-trace part (SSSS) and the purely double-trace part 
{DDDD) grow as A^ 7 . This is in contradistinction to the free field case, where (SSSS) ~ N® 
while {DDDD) ~ N 8 (see eq. (3.9)). Therefore, at one loop the contribution of the double- 
trace operator appears to be more suppressed than in the free case, if compared to that of 
the single trace. The double-trace contribution can compete with the single-trace one if the 
mixing parameter is of order one. This is however ruled out by the free field considerations 
above which set a to scale faster than 1/yN. 

In conclusion, we see that if one requires the large N free field theory amplitude to 
coincide with the one for single-trace operators (and therefore with the constant part of the 
supergravity induced amplitude), then the mixing produces no effect on the leading large 
N one-loop amplitude and, as a consequence, on its degeneracy. Such a behavior follows 
the supergravity pattern: The redefinitions of the supergravity fields (corresponding to the 
field-theoretic mixing) produce no effect on the four-point amplitude, provided that it is of 
the regular (i.e., not of the extremal or sub-extremal) type [19]. 

4 The J\f = 2 insertion procedure and Feynman rules 

We would like to find the quantum corrections to the lowest component (at #1,2,3,4 = 0) 
of the four-point correlator (Q k ) and identify the corresponding conformal invariants T 
in eq. (2.10). Already in the one-loop case, but even more so at two loops, the most 
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efficient technique, in our opinion, is the insertion procedure in N = 2 harmonic superspace 
[20, 17, 21, 22]. 

The basic ingredients of the M = 4 SYM theory in terms of J\f = 2 superfields are the 
hypermultiplet q + (x, 8 + , 9 + , u) and the N = 2 SYM fieid strength W(x, 6). The distinctive 
feature of the former is that it is a Grassmann analytic superfield (i.e., it depends only on 
the harmonic U(l) projections 9 + ,8 + ) while the latter is a chiral superfield. 

One of the advantages of the J\f = 2 formulation of the theory compared to J\f = 1 is 
that the HM composite operators like Tr(g + ) fc need no covariantization (no presence of the 
gauge superfield at the external points of the amplitude). Further, the HM matter interacts 
with the gauge sector only through a single cubic vertex. The true non-Abelian nature of 
the theory is encoded in the gauge self-interactions (as well as in the ghost sector, but we 
do not see it at the level <? 4 ). 

The quantum corrections to the four-point correlator (Q k ) at order g 4 can be obtained 
by a double insertion of the N = 2 SYM action 

Sm=2 sym = J d 4 xd 4 C , C = ^2 W 2 . (4.1) 

More precisely, 

(Q {k) ) 9 * = \f d 4 x 5 d 4 # 5 d 4 x 6 d 4 6 <Q«|£(5)|£(6)) tree . (4.2) 

Here (Q^\C(b)\C(Q))t rcc is a new, six-point correlator calculated at "tree level" (Born 
approximation) . The corresponding Feynman graphs are obtained by drawing the usual two- 
loop (order c/ 4 ) graphs for (Q^) and then inserting the linearized M = 2 SYM Lagrangian 
C into the gluon lines twice, each time into a different gluon line. 7 

The most important property of the new six-point amplitude is that it can be written 
in the factorized form 

(Q {k) \C\C) = G x A^ k - 2 \x,u) , (4.3) 

where Q is a particular nilpotent six-point superconformal covariant carrying harmonic U(l) 
charge 2 at points 1 to 4, and R charge 2 at points 5 and 6 (the chiral field strength W has 
R charge 1 while the HMs carry harmonic U(l) but no R charge). Since the left-handed 
9 a are the only superspace coordinates with positive R charge 1/2, we conclude that the 
expansion of G must start with four such 9s. The rest of the six-point correlator (4.3) is 
given by the function A^ k ~ 2 \x,u) which carries U(l) charge k — 2 at points 1 to 4 and can 
be expanded in the propagator basis as follows: 

k-2 

A( k ~ 2 \x,u) = X m Y k - m - 2 A m {x) . (4.4) 

m=0 

7 This insertion procedure is based on the formula for the second derivative of the four-point cor- 
relator with respect to the gauge coupling constant: (d/dg 2 ) (Q^) = / d 4 xsd 4 6s (Q^\C) + 
\ J d 4 £5d 4 6>5 d 4 :r6d 4 #6 (Q^ k ^\£\£). If both insertions are made into the same gluon line, this means in- 
serting the chiral-to-chiral propagator (W / (5)W(6)} ~ 5(5,6) into that gluon line. By performing the chiral 
superspace integration over point 6, one can show that the five-point correlator arising in this way precisely 
cancels against the single-insertion term in this formula. 
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The coefficients A m (x) are conformally covariant functions of the six space-time points. 
They have vanishing U(l) charge and are thus harmonic independent. 

The structure of the nilpotent covariant is determined by superconformal symmetry 
combined with the Grassmann analytic (or ^-BPS) nature of the first four points and the 
chiral nature of the last two points. The explicit form of (for 0+ = 0, r = 1, . . . , 4) has 
been worked out in Ref. [22]: 

tt4 2 2 2 2 4 4 
q _ llr=l x r5 X r6 x 12 x 34 x 13 x 24 ^ 

[12] 2 [34] 2 ri 4 r 23 + [14] 2 [23] 2 r 12 r 34 + [12] [23] [34] [41] [ti 3 t 24 - r 12 r 3 4 - ti 4 t 23 ] } . 
Here 



A{p r p s )(o r o s ) + p 2 o- 2 + p 2 s a 2 (4.6) 



and 



p r = (0+ - eiu+^i , a r = (0+ - eiu+)x^ , r = 1, . . . , 4 (4.7) 

are Q-supersymmetry invariant combinations of the analytic 9^ with the chiral 8^ 6 pro- 
jected with the SU(2) harmonics {u r )f . The harmonic contractions [12], [34], etc. and the 
polynomial Rj^=2 have been defined in eqs. (2.9) and (2.11), respectively. 

The aim of our two-loop calculation is to determine the factor A^ k ~ 2 \x,u) in the six- 
point covariant (4.3) and then to substitute everything in the insertion formula (4.2). Since 
we are only interested in the lowest component of the four-point correlator (Q k ), we can set 
all the external 9s to zero, 9^T = = 0, r = 1, . . . , 4. In this case is rather simple: 

2 2 4 4 

@ e+=0 = 9l9l R M=2 . (4.8) 



x 



Consequently, the six-point correlator becomes 

(Q^\c\c) e+=0 = eiei x ^AfA Rm=2 A ik-2) M . (4 . 9) 



56 



Further, substituting this into the double-insertion formula (4.2) and performing the trivial 
chiral integrations over 65 q, we obtain the two- loop correlator 



<Q( fc V =0 = x\ 2 x\,x\A RM=2 I A^\x,u) . (4.10) 

J X 56 



Rewriting the amplitude in the form (2.10), we can read off the following expression for the 
conformally invariant coefficient functions J-^ at two loops: 

•^m,fe-m-2,o( S ' = X 12 X 34 X 13 X 24 / 4 A m (x) . (4-11) 

J x 56 

Now, the practical question is how to compute A^ k ~ 2 \x, u) (and hence A m (x) and 
J-(s, t)) from the corresponding set of two- loop Feynman diagrams. It turns out that instead 
of setting 9 + = 0, as required in the final expression (4.10), it is much more convenient to 
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do the computations with #5^ = 0. The knowledge of the structure of G (4.5) allows us to 
easily switch from one of these superconformal frames to the other. 8 In the new frame G 
becomes even simpler, 



Q e5 , 6=0 =n^ + ) 2 > ( 4 - 12 ) 

r=l 

so that eq. (4.9) is replaced by 

(Q( fc )|£|£), 5 , 6=0 = f[(9+) 2 A^){x,u) . (4.13) 

r=l 

Then it is clear that in working out the expressions for the various Feynman graphs we can 
concentrate only on the terms with the maximal number of external 9s. In particular, at 
order g 4 this choice removes all graphs which contain non-Abelian interaction vertices. For 
example, the Y-shaped gluon subgraph in Fig. 5a vanishes because it has two chiral ends 
at the insertion points 5 and 6 and one analytic end (the gluon without insertion); after 
setting #5 : 6 = we are left with too few left-handed 9s at the analytic gluon end to supply 
the required R charge 2. Similarly, the TTT block in Fig. 5a has three chiral ends (in 
fact, only two, points 6 and 7 should be identified) and two analytic ends; once again, the 
analytic 9s cannot provide the required R charge 3. 



1 



1 



5a: 



5* *6 



• • • 

5 6 7 



TTT 



5b: ^ ^ 4 

5 5 6 

T TT 

Figure 5. Building blocks of the Feynman graphs. 

As a result of all these simplifications our task is reduced to listing all tree level Feyn- 
man graphs made out of the two building blocks T and TT in Fig. 5b. Although these 
blocks contain interaction vertices and hence integrals, the latter are easily computed using 



By "superconformal frames" we mean that Q can be cast into one of these forms by means of a finite 
superconformal transformations. 
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the harmonic superspace Feynman rules and they produce very simple rational space-time 
functions [20, 22]: 



= (q + (l)W(5)q+(2)) = ^\ [2l~]p{ 

}> = "'(2vr) 6 ^ 12 

where p,a have been defined in (4.7) and [12~] = (ui)fe^(u2)j (cf. (2.9)). 

Notice the characteristic presence of negative-charged harmonics in both expressions 
(4.14) and (4.15). This has to do with the important issue of harmonic analyticity [15]. 
The point is that the free HM satisfies the massless field equation 

D ++ q + (x,9 + ,9 + ,u) = 0, (4.16) 

where 

D ++ = u +i ^— - 4i9 + a^9 + ^- (4.17) 

is the supersymmetrized harmonic derivative (the raising operator of the group SU(2) re- 
alized on the charges ± of the harmonics). This equation can also be viewed as a Cauchy- 
Riemann condition on the harmonic coset SU(2)/U(1), hence the name "harmonic analyt- 
icity". Yet another interpretation of eq. (4.16) is that it defines q + as the highest weight 
state of an SU(2) irrep of charge 1 (isodoublet). When the HMs interact with the gauge 
sector, the harmonic derivative in (4.16) is modified by a gauge connection. Still, the gauge 
invariant composite operators like Tr(g + ) fc satisfy the same equation with a flat harmonic 
derivative. So, such operators correspond to the highest weight state of an SU(2) irrep of 
charge k. 

As the simplest example of harmonic analyticity in action, consider the HM two-point 
function (propagator). It has to satisfy the Green's function equation 

D++(q^(l)q+{2)) =5(1,2), (4.18) 

so here harmonic analyticity holds up to contact terms. The solution to this equation is 
(for simplicity we set 62 = 0) 

<? + U)9 + (2)> = ~ ^ ^ • (4-19) 

(x 12 + 4i [ ^9fo-6+ 

This can be checked by expanding the denominator in O^Of , using the relation □i(l/xf 2 ) ~ 
5(^12) and the naive harmonic differentiation rules 

D++6+ = 0, £>^ + [l~2] = [12], D+ + [12] =0, (4.20) 

as well as the rule 

= 5(U1,U2) (421) 

which follows from the relation d/dz(l/z) ~ S(z). In the context of our four-point investi- 
gations we always set all 9 = 0. Therefore here we never encounter harmonic singularities 



4c/ fabcfcde r-._ _-| 2 



TT 1562 = (q + (l)W(5)W(6)q + (2)) = - ^ 2 (4-15) 
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of the type (4.21) and can safely apply the naive rules (4.20). Then harmonic analyticity 
simply means polynomial dependence in u + and no dependence in u~. This is the case of 
the propagator (4.19) which now becomes 

(g + (l)g + (2)) e - =0 = M. (4.22) 

The same must hold for any n-point correlator of ^-BPS operators: They must be polyno- 
mials in the harmonic variables u + at each point. 

Clearly, the expressions for the building blocks T (4.14) and TT (4.15) are not harmonic 
analytic because of the presence of 1~ and 2~. This, however, is not a problem: The various 
building blocks or even complete Feynman graphs are not expected to be harmonic analytic, 
much like they are not conformal covariants. It is only the sum of all graphs that has these 
properties (see examples in Appendix A). In what follows we shall often profit from the 
expected harmonic analyticity of the final result to greatly simplify our graph calculations. 



5 Two-loop four-point amplitudes 



In this section we work in the frame #5^ 
the form 



where the six-point nilpotent covariant takes 



{Q (k) W)e 5 



6=0 



r=l 



k-2 

£ X m Y k ~ 

m=0 



*A m (x) 



(5.1) 



Here the coefficients A m (x) are not all independent, as the crossing symmetry 2 <-> 3 of eq. 
(2.10) implies that A m {x) is equal to A^_ m _2( x ) with X2 and x% interchanged. In particular, 

for k even, Ak-2 is invariant under this crossing transformation. By using eq. (4.11) one 

2 

can see that these facts are in accord with the corresponding transformation property of 
•^m i-m-2 o( s ' Therefore, we need compute only (^rp + 1) coefficient functions if k is 
even and if k is odd. 

At two loops the correlator (5.1) is given by the sum of a number of Feynman graphs 
containing four interaction vertices of HM with gluons and two C insertions. Let us discuss 
the general structure of these diagrams. According to the rules discussed in section 4, they 
are made out of the building blocks T and TT. The same argument based on 9 counting 
allows us to discard other sets of vanishing diagrams. Take, for instance, the graphs involving 
a HM loop with the insertion of two gluon lines connecting the HM propagators. Such a 
loop produces the expression (TT) 2 ~ p 4 a 4 = (see (4.15)). For the same reason any 
graph containing a corner with all the incoming or outgoing HM lines being free vanishes, 
since this implies concentrating too many 9s at some other corner. Finally, disconnected 
graphs containing subgraphs in the form of the one- or two-loop corrections to the two-point 
function of protected 1/2 BPS operators vanish as well. 

With these rules in mind and observing that at order g the maximal number of HM lines 
involved in interactions with gluons is four, we can now draw all possible configurations. 
The planar interaction topologies relevant to the calculation of the functions A m (x) are the 
same as the ones shown in Fig. 12 for the k = 4 case (thin lines). 9 The complete graphs 

9 Beyond those of Fig. 12 there are also other interacting topologies which, however, play only an auxiliary 
role as it will become clear from our procedure of identifying the diagrammatic contributions to A m (x). 
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are then obtained by combining an interaction topology from Fig. 12 (thin lines) with the 
matching free HM frame from Fig. 6a. Each resulting topology is labeled by a pair of 
integers (p, q) (0 < p < 2m, < q < 2{k — m — 2)) which fixes the particular structure of 
the free HM lines. This procedure automatically takes into account the possible crossing 
symmetry transformations of the end points. 




Figure 6. a) The free HM frame (p, q) associated with the function A m {x). b) The HM 
frame common to all graphs contributing to A m {x). 



Evaluating the explicit expressions for each topology amounts to multiplying the building 
blocks (4.14) and (4.15) and using a simple Fierz identity for the two-component spinors p 
and a. In this way the contributions of the individual graphs are always expressed in terms 
of the variables p 2 r , a 1 and r rs . 

Since the building blocks contain negatively charged harmonic variables, each individual 
graph involves some non-analytic harmonic dependence. The final sum (5.1) must however 
be harmonic analytic and moreover, the coefficients A m are harmonic independent. Thus, 
all the terms in the graphs with non-analytic dependence are in fact spurious and should 
cancel out. Due to the huge number of diagrams the evaluation of each diagram and the 
proof of the actual cancellation of non-analytic terms can be rather intricate (see Appendix 
A for a sample calculation). 

However, the knowledge that the final result has to be analytic allows us to skip the 
actual computation of all the diagrams and directly identify a minimal set of graphs relevant 
to restoring the correct harmonic analytic structure of the final sum. This can be achieved by 
a suitable procedure of identifying the harmonics. In Appendix B we describe this procedure 
in detail when applied to operators of weight k = 3,4. The corresponding minimal sets of 
relevant graphs for the k = 3 case are given in Fig. 11, whereas for k = 4 they are listed in 
Figs. 13 and 14. 

For generic k, the procedure of identifying the harmonics goes along the same lines 
as in Appendix B. In this way, we are led to select the minimal set of relevant diagrams 
which contribute to the function A m . They are drawn in Fig. 7 where the labels (p, q) are 
explicitly indicated for each graph. 
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a 



(m,k-]-m) 




a 2 

(m,k—3—tn) 



(m,k—2—m) 







r 




i 




a 3 


a 4 


(m+l,k-2- 












V 










(m,k—2—m) 


(m,k—2—m 



d 

(m,k—2—m) 



Figure 7. Relevant topologies of the interacting subgraphs of the planar graphs 
contributing to the function A m . Every subgraph has to be multiplied by the free 

contribution from Fig. 6a. 



We note that the topologies of interaction lines appearing in the diagrams of Fig. 7 are 
the same as the ones for k = 3, 4 cases (see Figures 11-14 in Appendix B). Let us explain 
why this happens by first looking at the case 1 < m < k — 3. We observe that the diagrams 
in Fig. 7 all contain the common free HM frame depicted in Fig. 6b. This frame provides 
the harmonic factor 

/ \ m-l / \ k-m-3 

xm _i yfe _ m _3 ^ f[ 12 ][43]J [13] [42] (5.2) 

which can be pulled out from all the graphs. It is easy to see that what remains forms the 
set of graphs contributing to the function A\(x) for the amplitude of weight 4 operators 
(see Fig. 14). Therefore, we can exploit the results for k = 4 to select the minimal set of 
interacting subdiagrams which survive harmonic identification. The particular cases m = 
and its crossing symmetry partner m = k — 2 can be treated in the same way. For m = 
the graphs and a.\ in Fig. 7 are absent and pulling out a common factor Y k ~ 4 we are led 
to the set of graphs for the function Aq of the k = 4 case (see Fig. 13). 

The important conclusion we reach at this point is that the two-loop calculation for 
arbitrary k is reduced to that for k = 4. The reason for this is quite clear: Since at two 
loops at most four HM lines are involved in the interactions, for a sufficiently large k the 
majority of the HM lines play the role of "spectators". 

Let us explain why the diagrams in Fig. 7 with the labels (p, q) indicated do indeed 
provide the contributions to A m . Recall that the graphs which contribute to A m are the 
ones from which we can eventually extract the propagator factor x m Y k ~ m ~ 2 . We have 
already pulled out a common factor x m ~ 1 Y k ~ m ~ 3 . Leaving aside for the moment the 
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graphs of the type 7a we see that all the other diagrams contain an additional factor 
XY ~ [12] [43] [13] [42]. Pulling it out we restore the required overall factor x m Y fc - m ~ 2 . 
What remains from the graph is a chargeless combination of harmonics which may contain 
constant harmonic independent terms, but also non-analytic terms. However, once we have 
pulled out the required overall factor, the rest has to be harmonic independent and can be 
computed by identifying all the harmonics as explained in Appendix B. 

We still have to take care of the graphs which contain interaction subgraphs of type 
7a. For instance, for the interaction topology 7ai besides the common frame (5.2) we can 
also pull out an extra factor [12] [43] ~ X. This leaves behind the free factor [13] 2 which 
is different from the expected [13] [42] ~ Y. The missing factor Y is restored only after 
summing up 7ai with all of its analyticity partners. The partners of the graph 7ai needed 
to restore analyticity, e.g., at points 1 and 3 are shown in Fig. 8. In Appendix A we analyze 
in detail how analyticity is achieved in this case and how such a long calculation can be 
efficiently bypassed by the procedure of identifying the harmonics. 




Figure 8. The graph 8a (topology 7ai) and its partners restoring analyticity at points 1 
and 3 (all the graphs have the same bottom part). The common free HM frame has been 

pulled out as the prefactor x m ~ 1 Y k ~ m ~ 3 . 



The reader may wonder why the graph 8d, which has the topology 7b, has been included 
in the set 8 where it clearly plays the auxiliary role of an analyticity partner of the graph 
8a. In fact, this graph plays a dual role: It becomes the principal graph contributing to the 
coefficient function A m +i(x). Indeed, pulling out the common factor x rn ~ l Y k ~ m ~ s and the 
extra factor X 2 contained in the the graph 8d, we find 



rig - [13][l-3"](p^ + piaf) 724 - [24][2"4-](p^ + fip\ 



(5.3) 



This expression is chargeless, so in it we can identify all the harmonics, after which we are 
left with the analytic term T13T24. We can say that the analytic part of eq. (5.3) contributes 
to the function A m +i, while the non-analytic part combines with the other graphs in Fig. 
8 into a harmonic analytic expression which gives a contribution to A m . 

Having identified the relevant graphs, we now proceed to computing the six-point cor- 
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relator. The Grassmann and space-time structures of the graphs in Fig. 7 are 
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ry£-* ry ^ ry ^ ry£-t 

x 12 x 34 x 13 x 24 



T14T23 - Ti3T 2 4 - T12T34 - T2z{p\(j\ + P%V\) ~ Tl4(P3°2 + P2 a 3 



.2,2 



„2~2n 



+ 7-13(^4^2 + P2 (j a) + 724^3 + P3°l) + r 12(p 4 tf 3 + P3V4) + T 3^(pl^2 + P2^1 

Collecting everything we arrive at the following complete harmonic analytic expression for 
the six-point correlator: 



(Q^\£\£) 



fc-2 



ry 2 ry 2 /y>2 rva 2 

x 12 x 34 x 13 x 24 m=0 



(5.4) 



/-id 



T12T34 



+ 
+ 
+ 



C^nMal + ploi) 



|C^T14T23 + {^sC b rn - C^Ti 3 T24 + 
(C£ " ^)T23(P?^ 4 2 + p\o\) + (C£ - 

(C£ - C«)7i 3 (p|ai + p|a 4 2 ) + (Ct - 

(Ci - C%)tM4 + P&D + (Cm - ^)r 34 ( /0 2 a 2 2 + pla 2 )} . 

where C m denote the combinatorial factors associated to the various graphs. 

The expression (5.4) has been derived for 1 < m < k — 3. As was already mentioned 
above, the cases m = (and its crossing partner m = k — 2) are a bit special, as for m = 
the graphs a3 and a4 are absent, while for m = k — 2 we do not have ai and a2- Otherwise, 
the derivation of the relevant contributions for these cases follows the same steps as above. 
Thus, the formula (5.4) remains valid for m = and for m = k — 2 as well, provided we 
impose the following formal requirement on the combinatorial coefficients 



Co 3 



°0 



°fc-2 



°fe-2 



0. 
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Before we proceed, let us briefly discuss the general symmetry properties of the combi- 
natorial factors. We note that for the graphs ai and &2 the range of indices is < m < k — 3, 
while for a3 and it is 1 < m < k — 2. For a generic m the combinatorial coefficients C m x 
and C^ 2 corresponding to the graphs a.\ and a2 are equal, C m x = C™ 2 = (though their 
analytic structures are different). Analogously, C^ = C^ 4 = C m and C^ = = C m . We 
also have the following additional relations 

Cm = Ck-m-2 > = Cfc-m-2 > (5-5) 

which in the special case m = m_ reduce to 

C m _ = ) Cm_ = ^m- ■ (5-6) 

Finally, to recast the six-point amplitude (5.4) in the form (5.1) we can use the relations 

Q+2 o+2 2 r 2 

2 2 _ _ fl+2,0+2 x rs x 56 7 N 

^ ~~ ~2 ' a r - 2 ' Trs ~ ^ ^ „2 „.2 „2 ™2 U 

valid for 65^ = 0. The very last step is to perform the space-time integration over the 
insertion points 5 and 6 in eq. (4.11). This introduces the functions $>W and $( 2 ) defined 
in (2.19) and (2.20) and their point permutations. Note that is totally symmetric under 
permutations of the external points. 

With these definitions at hand, from eqs. (5.1) and (4.11) we can identify 



<L m - 2 ,oM) = + (2C* - Ci)s + (2C m - CD) [*«(M)] 2 (5.8) 

+ (Ct - ci)UW{t/s, l/s) + (C d m - C%)*®(8,t) + (C m - Ci)-^ 2 \s/t, 1/t) 
s t 

Here x is an overall normalization coefficient independent of m, which is related to the 
normalization of the ^-BPS operators. 

The final result (5.8) essentially depends on the values of the combinatorial coefficients 
C m - The particular crossing symmetry properties (2.17) of the functions T± provide some 
additional relations among the corresponding combinatorial coefficients. Indeed, for m = 
m + we find 

C m+ = 2C m+ = 2C m+ (5.9) 



but there are no new restrictions on C m _ . 

We remark that the case k = 2 is special as there are no graphs of type 6a and the 

— O — ZUq. 



higher crossing symmetry requires 2Cq = Cq = 2C$. 
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The general expressions for the combinatorial coefficients C m valid for finite ./V are 



.a k -iP)(ai...a m b m+ i...bk-3qrn) 



amb m +i—bk-2Qn) 



O'mbrn+l ■ ■ ■ 



where 

(ai ...a k )= Tr(i (oi . ..t 0fc )) 

is the symmetrized (without l/k\) trace of fc generators of the color group. In these formulae 
the combinatorial factors are necessary to prevent the overcounting of the HM lines induced 
by symmetrization. 

To compute these coefficients in the large N limit we use the conventions [t p , t q ] = ifp qe te 
and fabefabe' = 2N5 ee t as well as the following fusion and splitting rules 

(a 1 ...a l A)(a 1 ...a l B) « (I + l) 2 /!iV'" 1 Tr(A5) , 
Tr[ai...aia 1 ...aiA] « N l TrA. 

In the large N limit we find 

C^ = 2C b m = Ci = k 4 N 2k , C c m = 2k A N 2k . (5.10) 

In this formula the range of the index m is < m < k - 3 for C« , < m < k - 2 for , 
1 < m < k - 2 for and < m < k - 2 for C^. The case of Cg deserves some special 
attention. Looking at the graph 6b for m = 0, we realize that the gluon lines can "slide" 
along the HM loops, thus giving rise to two equivalent planar diagrams. Hence, the factor 
Cq is twice as big as the other with m/0: 

C h = k*N 2k . (5.11) 

So, the combinatorial coefficients have the remarkable property that they do not depend on 
m, except for the single "anomalous" coefficient which changes its value only for m = 0. 

Substituting these combinatorial coefficients into eq. (5.8), choosing canonical normal- 
ization for the ^-BPS operators and recalling the definition (2.16), we obtain our main 
result (2.18). 
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6 Two-loop four-point amplitudes from the OPE 



Here we show that the four-point two-loop amplitude for ^-BPS operators of weight 3 can 
be completely reconstructed just on the basis of the expected operator product expansion 
structure and the knowledge of the one-loop anomalous dimensions for certain operators of 
twist 2 and 4. Analogously, we obtain some restrictions on the structure of the amplitude 
for weight 4 operators. The discussion in this section provides an independent check on our 
diagrammatic computation. 

6.1 Weight 3 

The diagrammatic treatment clearly shows that at two loops the four-point amplitude is 
made out of two basic conformal integrals, (2.19) and (2.20). Therefore, the most general 
Ansatz compatible with the symmetry requirement ,F(s,t) = 1/t T(s/t, 1/t) is 

^M) = ^•^(es + riit + i))^^)] 2 (6.1) 



4 



+ ^ 2 \t/s,l/s) + u(^ 2 \s,t) + -$ {2) (s/t,l/t) 

s \ t 

i.e. each function depends on four unknown coefficients g, r], fi and v. 

Already at the level of the free field theory one can see that the OPE O^O^ of two 
tj-BPS operators exhibits a heredity property: If some superconformal primary operator 
emerges in the spectrum of this OPE then it must also appear in the spectrum of OWoW 
with n > k. Therefore, knowing the spectrum of the operator product 0( 2 )C?( 2 ) (see Refs. 
[25, 23]), we can use it to obtain restrictions on the four-point amplitude underlying the 
OPE of the ^-BPS operators with higher weights. From the technical point of view, we 
need to build the conformal partial wave amplitude (CPWA) expansion of the conformal 
integrals entering eq. (6.1) and then to match the corresponding short-distance expansion 
of the four-point amplitude with perturbative CPWA contributions of operators with known 
anomalous dimensions. In this way one can fix the undetermined numerical coefficients in 
eq. (6.1). 

The CPWA expansions of the conformal integrals in eq. (6.1), corresponding to taking 
a short-distance limit x\ — > X2, x$ — ► x±, have already been obtained in Ref. [23] (cf. 
Section 3 there). They are formulated in terms of the conformal cross-ratios v = s/t and 
Y = (t — l)/t. Under the assumption that the conformal dimension A of an operator is 
of the form A = Ao + 7, where Ao and 7 are the canonical and the anomalous parts, 
respectively, the non-analytic term v^^ 2 entering the CPWA of this operator gives rise to 
perturbative logarithms 

v A/2 = i + i/ 2 7 lnu + l/8 7 2 ln 2 u + ... 

Therefore, the information about the one-loop (order A) anomalous dimensions is reflected 
in the coefficients of the In 2 v terms in the perturbative CPWA expansions and this is how we 
can use it to derive the form of the two-loop (order A 2 ) amplitude. The explicit calculations 
are not particularly instructive and therefore in the following we mainly restrict ourselves 
to the formulation of the results obtained. 

Operators of twist 2. 
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The R-symmetry sing let of the OPE O^O^ contains an infinite tower of twist 2 oper- 
ators of increasing spin, the lowest member being the Konishi scalar of canonical dimension 
2 and of one- loop anomalous dimension 3A. The requirement of reproducing this field in 
the OPE derived from eq. (6.1) fixes two of the four unknown coefficients, namely, 

3 2 

ri = v = — . (6.2) 

Given relation (6.2), the whole tower of twist 2 operators is then correctly reproduced by 
the four-point amplitude. However, the operators of twist 2 produce no restrictions on the 
coefficients g and v. 

Operators of twist 4 in the singlet channel. 

To get further restrictions on the four-point amplitude we consider the superconformal 
primary operators of twist 4 in the singlet R-symmetry channel. The simplest example are 
the four quadrilinear dimension 4 operators £12, £± studied in Ref. [24] (see also [5]). These 
operators diagonalize the matrix of the dilatation operator at one loop with the following 
result for the large TV one-loop anomalous dimensions 10 

71 = 0, 7 2 = 6A, 7£ ± = i(13±>/4l)A. (6.3) 

By using the explicit form of the (canonically normalized) operators £12, £± established in 
Ref. [24] we then compute their free three-point functions with the canonically normalized 
Tj-BPS operator in the large N limit and find 

while (0( 3 )C>( 3 )£ 

1,2) ~ 0(1/N 2 ). Such a suppressed behavior of the three-point amplitudes 
involving £12 is naturally explained by the fact that these fields are double-trace opera- 
tors. Thus, among the four operators only £± participate in the large ./V operator product 
expansion of the four-point amplitude. The normalizations of their CPWAs are given by 
the constants A±. 

Combining the free, the one-loop amplitudes [7] and the two-loop Ansatz (6.1), we have 
worked out the corresponding CPWA expansion and have found the following system of 
equations 11 encoding the large iV contribution of the fields £±: 

A + + A_ = , (6.5) 
^ +7S+ +A_ 7S _ = §^2> (6-6) 

a 2 , 2 / 81 8 4 \ ^ fo 

A ^ +A -^- = {^ + i e+ -^)w- (6 ' 7) 



10 At fixed A the anomalous dimension 71 vanishes only when N — 00 [24]. 

11 These equations are obtained after careful disentangling of the contribution of the superconformal de- 



scendants of the twist 2 fields. 
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It is easy to check that the first two equations are indeed satisfied with the values of the 
A's and 7's found, while the third equation yields 



Q = -\- (6.8) 

At this stage only one coefficient, e.g. /x, remains undetermined. 
Operators of twist 4 in the irrep [0,2,0]. 

The OPE CK 2 )0( 2 ) has ten different SO(6) channels. Among them only three may 
contain unprotected superconformal primary operators: [0,0,0], [0,2,0] and [1,0,1]. As we 
have already exploited the known information about the lowest dimensional operators in 
the singlet, now we move to the superconformal primaries of dimension 4 arising 12 in the 
irrep [0,2,0]. 

There are four superconformal primary operators of canonical dimension 4 which might 
contribute to the operator product expansion under study. One of them is the protected 
semi-short double-trace operator [25]. The other three are long unprotected operators O, Q± 
whose one- loop anomalous dimensions were found in Ref. [26] (see also [27]): 

70 = 3A, 7 e ± = ^(5±VE)\. (6.9) 

The protected operator as well as O are double-traces, therefore we do not expect them to 
contribute to the four-point amplitude in the large N- limit. Denoting by B± the normal- 
izations of the CPWAs of the operators 0± and proceeding in the same manner as before, 
we find the following system of equations: 

147 

B + +B - = m^ (6 - 10) 

B +7e+ +B-7e- = ^, (6.11) 
35A 2 

B +1 2 e+ + B^ 2 e _ = —(9 + 2fi) . (6.12) 
The first two equations are used to find 

21 

£± = ^(7^375). (6.13) 

The normalization constants are positive, as they should (the squares of the normalization 
coefficients of the three-point functions). Substituting these values in the third equation we 
find fi = 0. 

Thus, by exploiting the information about the one- loop anomalous dimensions of some 
superconformal primaries, we have completely determined the freedom in the most general 
Ansatz (6.1): 

32 

g = H = 0, r} = v = —. (6.14) 



12 There is only one ^-BPS operator of dimension 2 in the irrep [0, 2, 0] and it does not produce any new 
restrictions on the four-point amplitude. 
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One can easily see that with these values the amplitude (6.1) indeed coincides with the 

(3) 

function !F\ found in the previous section. 

Finally, we remark that in Ref. [5] the two-loop anomalous dimensions for the operators 
0± were obtained 

17±5x/5 2 

7± = - A . (6.15) 

Working out the logarithmic terms in the two-loop OPE implied by the four-point amplitude 
for k = 3 we have checked that eq. (6.15) is indeed compatible with our findings provided 
the one-loop corrections to the normalization constants are 

91 r 1 i 

B± = (0®0®Q±) 2 = ^(7 T 3^5) [l - -(35 ± 3^5)a] . (6.16) 

6.2 Weight 4 

Now we obtain some restrictions on the four-point two-loop amplitude for the |-BPS op- 
erators of weight 4 from the knowledge of the one-loop anomalous dimensions for the same 
operators of twist 2 and 4 that were already discussed in the previous section. 
For k = 4 we expect two functions J-±^ with the symmetry property 

T^(s,t) = l/tT^(s/t,l/t). 
The general Ansatze compatible with these symmetries are 

4 4) M) = ^■[\(g±s + v±(t + i))[^ 1 Ks,t)] 2 (6.17) 

+ ^$ (2) (t/s, 1/s) + v ± (V 2 ) ( 8 t ) + 1$( 2 ) (s/t, 1/t) 
s \ t 

i.e. each function depends on four unknown coefficients g, rj, \x and v. 
Operators of twist 2. 

This time the Konishi field (as well as its higher-spin cousins) imposes the following 
relations 

4 2 

v+ = u + = -=A, (6.18) 

while the other coefficients remain invisible. Clearly, these values of rj+, v+ coincide with 
those in eq. (2.18) obtained by explicit diagrammatic computations. 

Operators of twist 4 in the singlet channel. 

A calculation similar to the case k = 3 produces the following system (after disentangling 
the contributions of the superconformal descendants of the twist 2 fields): 

A + + A_ = » , 

A h7 £ + +A_ 7S _ = ^A ) (6.19) 

, n n 2 / 6 \ A 2 

^+7s + + ^-7e_ = 3 (3^+ -- + ??_ + 2^. + 2v- + g-)j^, 
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where 7s ± are given by eqs. (6.3). Solving this system we find 



rj- = 12 - 3g+ - - 2v_ - g_ . (6.20) 
This completes our consideration of the singlet channel. 
Operators of twist 4 in the irrep [0,2,0]. 

The occurrence of the superconformal primaries Q± in the double operator product 
expansion implied by the four-point amplitude leads to the following system of equations 

Bi+B _ 2744 



25N 2 ' 

784A 

B +le+ + 5_ 70 _ = — T , (6.21) 



5N 2 
98A 2 

5+7© + + B-7e_ = (r/- + 2^_ + + q. 



It allows us to find 



ry_ = 12 - 2/i_ - 2^_ - g_ . (6.22) 

Comparing the last formula to eq. (6.20) gives g + = 0, which coincides with the corre- 
sponding value of the coefficient g + in eq. (2.18). 

This exhausts the predictive power of the method. The remaining coefficients fi±, V- 
and p- can be found only from the diagrammatic computation. If we had the information 
about some other superconformal primaries appearing in the OPE for k = 3 or k = 4, we 
could use it to further refine our Ansatz. 

Finally, according to formula (2.18) we have 

fi + = ij- = v- = , /x_ = £>_ = 4 . (6.23) 
With these values eq. (6.22) is trivially satisfied. 
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Appendices 



A Harmonic analyticity constraints 

In this Appendix we develop a set of systematic rules to select the relevant graphs con- 
tributing to a correlation function and find consistency conditions for the combinatorial 
coefficients. These rules are based on the general property of harmonic analyticity satisfied 
by any correlation function of composite operators in G-analytic harmonic superspace. 

Harmonic analyticity for a composite operator O k = Tv(q n q k ~ n ) means D ++ O k = 0, 
where D ++ is the harmonic derivative given in (4.17). Correlation functions then satisfy 

D++(O k (l)---O k (n)) =0 i = l,---,n (A.l) 

(modulo contact terms). At every order in perturbation theory we can implement this 
condition directly on Feynman graphs. As a consequence of the identity 

D++(<t(l)q + (2)) = -{q+ (l)q + (2))D+ + = 5(1,2) (A.2) 

where 5(1,2) is the delta function in the G-analytic harmonic superspace, every time the 
derivative hits a HM propagator shrinks it to a point. In particular, when D ++ hits a free 
line we have a vanishing contribution (we neglect contact terms), whereas when it hits a 
HM line involved in an interaction with a gluon, the effect is to shrink the propagator line 
to a point and move the internal gluon vertex to coincide with an external vertex. A simple 
example is given in Fig. 9 where analyticity implemented at the vertex 1 of diagram ai 
produces the shrunk configuration a2- 

If the same shrunk configuration can be obtained when hitting different graphs (but at 
the same point), we can say that all these graphs conspire together to restore harmonic 
analyticity at that point. Repeating this procedure, one eventually recovers the whole class 
of graphs needed to form a harmonic analytic amplitude. 



Figure 9. Example of an analyticity condition. The arrows point at the propagator hit by 



D++. 



Again, we exemplify this by taking the simple case of Fig. 9. There, the same shrunk 
configuration is obtained both from diagrams ai and a3. Since these are the only two dia- 
grams which give the configuration &2 when hit at the point 1, the sum of their contributions 
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must result into an expression which is analytic in the 1 variable. In fact, using eq. (4.14) 
their contributions read 

ai - C ai [13][42] {[21-}pl + \\2~\p\ - 2( Pl p 2 )} {[43"]p| + \U~\pl - 2{ PzPi )} 
a 3 - C 03 [12][42] { [31-]p? + [13"]^ - 2( Pl p 3 )} {[43"]pl + [34"]^ - 2{ Pm )} 

(A.3) 

where C ai are the combinatorial factors (with possible signs included). The two expressions 
are non-analytic in the 1 variable because of their nontrivial dependence on . The action 
of Df + on the sum of the two contributions must give zero, according to the condition 
(A.l). Remembering that D± + only acts on the u± harmonic converting it into uf , we 
obtain 

(C ai +C 03 )[42][13][21]{[43-]p?pi + [34-]p?pl} = (A.4) 

which implies 

Cai = ~Ca 3 (A-5) 

Using this condition, the sum of the two contributions in (A.3) reads (we consider only the 
part depending on 1 _ ) 

[42] {[13] [21-] - [12][31-]} {[43-]p?p| + [34"]^} (A . 6) 
The harmonic cyclic identity 

[13][21"]-[12][31-] = [23] (A.7) 

shows that the non-analytic contributions from the two diagrams cancel against each other 
leaving an analytical result in the 1 variable. 

As appears from this simple example, imposing analyticity conditions allows to draw 
systematically all the diagrams which at a given order contribute to the analytic final result 
for the correlation function and to find nontrivial relations like (A. 5) which constrain the 
combinatorial coefficients. These relations can be used as a check when the coefficients are 
directly computed from the diagrams. 

We notice that the final analytic result can be obtained by using a shortcut procedure: 
Given the sum of the two 1 _ dependent contributions in (A.3) we make the harmonics 
identification 1 = 2. As a consequence, the contribution from a 3 vanishes whereas ai gives 

C B1 [42][13][21-] {[43-]p?p! + [34"]^} =H C 0l [42] [23] { [43"]^ + [34"]/^} (A . 8) 

which is the correct result. The identification of the harmonics we chose is not the only 
possible one. We could have made the alternative identification 1 = 3. In this case the 
contribution from diagram ai vanishes whereas from diagram a 3 we obtain 

C as [12][42][31-] {[43-]p?p| + [U-] P \ P l} ^ -C as [23] [42] {[43-] P \ P \ + [34"]/^} (A . 9) 
which again is the correct result. 
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We now consider the cases of our interest. We start by studying the analyticity condi- 
tions for the set of diagrams of Fig. 8 in the main text. They have in common two free 
lines which give the analytic structure [12] [43] and a "base" T245T246: 

T 2 45T 24 6 = [24-][42-](p|ff| + p\a\) + 4(p 2 p 4 ) (<t 2 <7 4 ) 

= r 2A -[24][2-A-}(pla 2 2 +p 2 2 al). (A.10) 

It is easy to see that leaving aside the base as well as the free propagators leads to the 
following harmonic and Grassmann structures: 



(8a) 
(86) 
(8c) 

m 



[13] 2 [T 125 T 4 36 + 5 ^6] = [13] 2 [21-][43-](/>?a 3 2 + plaj) 
[13][43][T 125 Ti36 + 5^6] = [13][43][21-][13-](p?<7| + p 2 3 al) 
[13] [12] [T 435 T 136 + 5^6] = [13][12][43-][31-](p? ( rl + p| ( 7?) 
[12][43][r 135 r 13 6 + 5^6]= 2[12][43] {-[l^ ^]{p\al + pfa 2 ) + r 13 } 

(A.ll) 



The factor 2 in the graph (8d) is due to the symmetrization 5 <-> 6. Note that this is the 
only graph from this set which gives an analytic contribution proportional to ri 3 . The other 
Grassmann term p\o\+p\(j\ always appears with non-analytic harmonic factors. Collecting 
all such terms we find 

[13] {C a [12] [43~] [31-] + C 6 [13][21-][43"] + C c [43][21-][13-] - 2C d [12][43][l-3"]} (A.12) 

where we have introduced combinatorial factors C a , ■ ■ ■ ,Cd- We want to achieve harmonic 
analyticity at points 1 and 3. According to the rule described above, we hit the graphs with 
harmonic derivatives Df + and Dt . 

We begin by applying Df + to the expression (A.12) to obtain 

[13][12]{(-C 6 -2C d )[43][13-] + (C o + C c )[43-][31]} = (A.13) 

Since the two harmonic structures within the brackets are linearly independent, we have 
the following two conditions 

C b = -2C d C a = -C c (A. 14) 

Similarly, hitting (A.12) with D~£ + we obtain 

[13] [43] {(C a + C 6 )[21-][13] + (C c + 2C d )[12][31-]} = (A.15) 

which implies 

C c = -2C d C a = -C b (A. 16) 

The last set of conditions is consistent with (A. 14) if C = C c . This was somehow expected 
since graphs 8b) and 8c) only differ by a reflection (1 <-> 3, 2 «-> 4) which should not affect 
the combinatorial factor. 

We can solve the constraints (A. 14) ending with only one arbitrary coefficient, e.g. C d . 
With this choice for the coefficients the expression (A.12) becomes harmonic analytic with 
the help of the identity 

[43][21-][13-] + [12] [43-] [31-] - [13] [21-] [43-] + [12] [43] [1" 3"] = [42] (A.17) 
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Finally, taking all this into account, attaching the "base" and restoring the missing 
space-time factors, we find the total contribution of the graphs in Fig. 8 to be 

Total Fig. 8: ff 2] 4 [43 j {-[13][42](p 2 a 3 2 + p 2 a 2 ) + [12][43]r 13 } 

x 12 x 13 x 24 x 34 

x {-[24][2-4-](^o| + plal) + r 24 } (A.18) 

We see that analyticity at points 1 and 3 has indeed been achieved but not at points 2 and 
4, of course. 

As in the previous example, the correct result can be obtained by using harmonics 
identification. In fact, pulling out the free line contribution [12] [43], we can identify the 
harmonics 1 = 2 and 3 = 4. As a consequence, in eq. (A. 12) the contributions from 
diagrams a,c,d vanish and we are left with 

a[13] 2 [21-][43-](p^3 2 + pla\) X m" 4 C b [13] 2 (pla 2 3 + p 2 a 2 ) . (A.19) 

Since the final result must be analytic, we can find an analytic term which under the above 
identification brings to the same expression 

-C 6 [13][42](p? o\ + p\a\) C 6 [13] 2 (p 2 <x 2 + p%a\) . (A.20) 

Therefore this must be the correct term which appears in the final result and it is indeed 
what we found in (A.18). 

We notice that using this shortcut we can easily figure out the complete analytic struc- 
ture for the sum of diagrams in fig. 8. In fact, analyticity at points 2 and 4, which should 
be restored by adding another set of partner graphs, can rather be obtained by identifying 
2 = 4. This identification leaves the single term T24 in the bottom part of eq. (A.18) which 
becomes completely analytic. 



We now consider the more complicated situation of graphs in Fig. 10. By imposing 
analyticity at points 1 and 4 we can find consistency conditions between the coefficients of 
the graphs 10c) which have been checked independently in Section 6, and the coefficients 
10a) and 10b) which don't have any independent check. These conditions will necessarily 
involve new graphs contributing only to non-analytic structures, necessary in order to cancel 
similar terms coming from the graphs 10a), 10b) and 10c). 

We start by considering the graph a) and apply the outlined procedure in order to 
achieve harmonic analyticity at points 1 and 4. It is easy to realize that in the large N 
limit, i.e. neglecting non-planar graphs, the set of conspiring graphs is just the one showed 
in Fig. 10. 

As in the previous example, we write the harmonic and Grassmann factors coming from 
these graphs: 
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f , f , f , 

Figure 10. Graphs conspiring to harmonic analyticity at points 1 and 4. 

(a) : [12] [13] [42] [43] [T 135 T 126 T A25 T A36 + 5 ~ 6] ^ 

[12] [13] [42] [43] { [[31-] [12-] [24-] [43-] + [13~] [21"] [42"] [34"]] 

r 2 2 2 2 i 2 2 2 21 1 

(b) : [12] [13] [42] [43] [T 136 T 125 T A25 T A36 + 5 ~ 6] ^ 

[12] [13] [42] [43] { [31-] [12-] [24~] [43~] [p\a\p\al + a\p\alp 2 A \ 

+ [13-] [21-] [42-] [34-] [g\gIp\ P 1 + p\pla 2 3 a 2 A ] } 
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(ci) 


[i2] 2 [43] 2 [r 1 3 5 r 1 3 6 r 425 T4 26 + 5^6]-^ 






2[12] 2 [43] 2 {-[13][l-3-](p 2 ( T| + p| ( T 2 ) + r 1 3}{-[24][2-4- 


](picr| + P4<T2) +t 24 } 


(c 2 ) 


same as (ci) with 2 <-> 3 or 1 <-> 4 




(di) 


[12] 2 [13][42][43] 2 [(Tr) 15 63(Tr)4652 + 5 ^6]^ 






[12] 2 [13][42][43] 2 [l-3-][4-2-]( /0 2 p 2 a 3 2 a 2 + a\a\p\p 2 A ) 




(4) 


same as (cfi) with 2 <-> 3 or 1 <-> 4 




(ei) 


[12][13][42][43] 2 [ri35Ti26(Tr)4562 + 5 ^6] — ► 






[12][13][42][43] 2 [2-4-][13-][21~](a 2 ( r 2 p 2 p 2 + p 2 /^ 2 ) 




(e 2 ) 


same as (ei) with 2 <-> 3 




fei) 


same as (ei) with 1 <-> 4 




(e 4 ) 


same as (ei) with 2 <-> 3 and 1 <-> 4 




(/i) 


[12][13][43] 2 [r 1 35r 126 r 42 5T4 26 + 5^6]-^ 






[12][13][43] 2 [13-][21-](a?p 2 + p 2 <7 3 2 ) {-[24][2-4-](p 2 a 4 2 + 


Pi "!) + T 24} 


(/ 2 ) 


same as (/i) with 2 <-> 3 




(A) 


same as (/i) with 1 <-> 4 




(/4) 


same as (/i) with 2 «-> 3 and 1 <-> 4 


(A.21 



We concentrate on the contributions to the three independent structures not containing r 
factors 



tt — 2222, 2222 
■tr — 2222, 2222 

Z = p\ololp\ + o\plp\ol (A.22) 

They are 

(a) : C a [12] [13] [42] [43] {[31-] [12-] [24~] [43"] + [13-][21"][42-][34-]}z 

(b) : C 6 [12][13][42][43]{[31-][12-][24-][43-]V+ [13-][21-][42~][34-]?7} 

(ci) : 2C c [12] 2 [43] 2 [13][l-3-][24][2-4-](f/ + Z) 
(c 2 ) : same as (ci) with 2 <-> 3 or 1 «-> 4 

(di) : C rf [12] 2 [13][42][43] 2 [l-3-][4-2-][/ 
(d 2 ) : same as (di) with 2 <-> 3 or 1 <-> 4 
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(ei) : C e [12][13][42][43] 2 [2-4-][13-][21-]C7 

(e2) : same as (ei) with 2 <-> 3 

(es) ■ same as (ei) with 1 <-> 4 

(e^) : same as (ei) with 2 <-> 3 and 1 <-> 4 

(/i) : C / [12][13][43] 2 [13-][21-][42][2-4^](t/ + Z) 

(/2) : same as (/1) with 2 <-> 3 

(fs) : same as with 1 <-> 4 

(/4) : same as with 2 3 and 1 <-> 4 (A. 23) 

We have assumed that diagrams which differ by the exchange 2 <-> 3 or / and 1 <-> 4 have the 
same combinatorial factor. 

By acting with Df + and setting to zero the independent structure proportional to U 
we get 

[12] [13] [42] [43] {[12] [43] [13-] [2-4-] ( - 2C C - C d - C e - C } ) 

+[12][34-][13-][42-](-C 6 -Ce-C/)} = (A.24) 

Doing the same for the structure proportional to V gives an equation like (A.24) with 2 <-> 3, 
then it will not impose any new condition on the coefficients. The part proportional to Z 
gives instead 

[12] [13] [42] [43] { [12] [43] [13"] [2~4T] ( - 2C C - Cf) 

+[12][34~][13-][42-](- C a - C/)} + 2 3 = (A.25) 

By now setting to zero the two linearly independent analytic structure appearing in (A.24) 
and (A.25) we get the following conditions relating the coefficients of the graphs of Fig. 10 

C a = 2C C = — Cf 

C b = C a -C e = C a + C d (A.26) 

We notice that the condition Cf = —2C C has been already obtained in (A. 16) (diagrams 
10c) and lOf) are the same as 8d) and 8c), respectively). 

It's now easy to see that there are no extra conditions coming from the terms containing 
one t factor. The same is true if we impose the analyticity condition at point 4. So we 
conclude that (A.26) is the full set of conditions imposed by the requirement of partial 
analyticity at points 1 and 4 on the graphs in Fig. 10. 

Using the relations (A.26) one can check that the total contribution obtained by sum- 
ming the expressions in (A. 21) ends up to be analytic in the 1 and 3 variables. 

Again, as in the previous examples, the final answer can be obtained by suitable iden- 
tification of the harmonics. 
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To summarize, what we have leaned from the examples described above is the following: 
On one hand, imposing harmonic analyticity on the final result allows to find consistency 
conditions among the coefficients of different sets of graphs. On the other hand, once we 
give analyticity for granted, we can bypass the whole procedure of evaluating all diagrams 
conspiring to a final analytic result by performing harmonics identification, so drastically 
reducing the number of graphs to be computed. 



B Sample calculations: k = 3 and k = 4 

In this Appendix we describe in detail the calculation of the A m functions (see eq. (5.1)) 
for the cases k = 3 and k = 4. In particular, we show how the procedure of identifying the 
harmonics can be used both to unambiguously select the diagrams which contribute to a 
given function and to drastically reduce the number of diagrams one needs compute. 



B.l The case k 



The expected form of the complete amplitude for k = 3 case is 

■[13] [42] 



(Q (3) |£|£> 



2 

x 13 x 42 



(B.l) 



As discussed in the main text, the two functions Aq,A\ are related by crossing symmetry, 
and we need compute only one of them. To calculate for example Aq, it is sufficient to 
look at contributions proportional to the harmonic structure [13] [42] in (B.l). They can be 
unambiguously selected by identifying the harmonics 1 = 2 and 3 = 4 in (B.l) 



{Q^\c\c) t 



'5,6 = 



[i3]; 

X 13 X 42 



-A (x) 



(B.2) 



When looking at perturbative contributions at two loops, the harmonics identification im- 
plies the vanishing of all graphs containing: 

(i) at least one free line 1 — > 2 or 3 <— 4 

(ii) the blocks TTi 562 or TT 356 4 

(iii) the blocks T125T125 and alike 

Discarding also the graphs which vanish by simple theta counting (see the general discussion 
in Sect. 5), the relevant surviving diagrams are then shown in Fig. 11. Upon the above 
harmonic identification, they already contain the analytic factor [13] 2 required in eq. (B.2). 
Since the coefficient function Aq is independent of the harmonics, once we have pulled [13] 2 
out, we can go a step further and identify the harmonics 1 = 3 in the rest of the expression. 

As a consequence, the contributions from graphs 11c and lid drop out. In the rest the 
harmonic dependence reduces to a common factor and, neglecting the space-time structure, 
the final contribution of each diagram is given by its combinatorics 

2Ci n) = 2 x 3 4 iV 6 



r (n) 



ci 11 ) = c 



(11) 
/ 



2C; 



(11) 



= 3 4 N 6 



2 x 3 4 N 



3 N 

6 



(B.3) 
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Figure 11. Relevant graphs for the k = 3 case. 

The extra factor 2 in front of ci 11 ^ and C^ 11 ^ is due to the symmetry of the corresponding 
diagrams under the exchange 5 «-> 6. 

In this simple example it is already clear how one can exploit the condition of harmonic 
analyticity of the final result. It is in fact this condition which allowed us to make the 
identification 1 = 3, thus reducing the number of Feynman diagrams to be computed. In 
a complete calculation performed without identifying the harmonics, the diagrams which 
vanished under this identification would contribute only to cancel non-analytic terms coming 
from the rest of the diagrams. 



38 



B.2 The case k = 4 



We now consider the more complicated case k = 4. There we have 



1 



<Q (4) |£|£k. 



(B.4) 



r=l 



X 




This correlator is crossing symmetric under the exchange 2 <-> 3 which relates the coem- 
cient functions and ^2 to each other, while A\ is crossing symmetric and independent. 
Therefore, we need to perform the calculation in two independent channels. 

To select the relevant graphs which contribute in the two channels we observe that in 
the k = 4 case we have 8 HM lines and 2 gluons, which can at most connect 4 HM lines. 
Therefore, at least 4 HM lines remain free. If these free lines form a free corner, i.e. they all 
come out from the same vertex, we know that such graphs vanish because of theta counting 
(see the general discussion in sect. 5) . It follows that the free lines have to form at least 
one disconnected pair. 

If the free pair is [13] [42] , such graphs can contribute to the functions Aq, A\ in eq. (B.4) 
but not to Ai- Then we can write down the contribution of all such graphs with a free pair 
[13] [42] in the form 



Pulling out the factor 1 ^ 2 ' is equivalent to removing the corresponding free lines from the 



graphs, which results in a configuration with k = 3. Then we can go a step further and 
identify the harmonics pairwise within the brackets. If we identify 1 = 2, 3 = 4 and the 
graph does not vanish, then it contributes to Aq. If instead we identify 1 = 3, 2 = 4 and 
the graph does not vanish, then it contributes to A\. Note that it may happen that the 
graph vanishes under each of these identifications, then it should be discarded. 
Alternatively, if the free pair is [12] [34], eq. (B.5) is replaced by 



Once again, removing the free lines we obtain a k = 3 configuration in which we identify 
the harmonics pairwise within the brackets. If we identify 1 = 2, 3 = 4 and the graph does 
not vanish, then it contributes to A\. If instead we identify 1 = 3, 2 = 4 and the graph 
does not vanish, then it contributes to A2. 

We now use this strategy to select the relevant diagrams which eventually contribute 
to the coefficient functions. We carefully draw all the diagrams which are nonvanishing for 
theta counting and select the ones which survive under one of the harmonics identifications 
described above. We are then led to the set of graphs on Fig. 12 where, in order to 
reduce the number of diagrams, we have indicated with thick lines all free HM propagators, 
whereas thin lines are interacting (they contain an interaction vertex and gluons connect 
these vertices in all possible ways). All the graphs have 4 free lines and they have been 
organized according to the number of free lines coming out of a single point. We have 




(B.5) 




(B.6) 



39 



not drawn those diagrams that vanish after harmonics identification (the analog of graphs 
11c and lid of the previous example). As already explained, in an exhaustive calculation 
done without identifying harmonics these diagrams would conspire to cancel nonanalytic 
contributions from the rest. 




Figure 12. Relevant graphs for k = 4. The thick lines are free HM propagators whereas 
thin lines contain vertices of interaction with gluons. 

We now apply the harmonics identification procedure to figure out which diagrams 
contribute to which function. Let us analyze in detail the graphs in Fig 12a. The first kind 
of graphs contain a free pair [13] [42], so they belong to the type of eq. (B.5). According 
to the rule above, we pull out the free factor and then we identify the harmonics pairwise. 
The identification 1 = 3, 2 = 4 annihilates the graphs, so they cannot contribute to A\. On 
the contrary, the graphs can survive the identification 1 = 2, 3 = 4, so they contribute to 
Aq. If we now draw the gluons, we see that only one configuration is allowed and the graph 
is reduced to lie from the case k = 3. The same argument applies to the second kind of 
graphs 12a, reducing them to a single structure like llf. Similarly, the third and the fourth 
graphs 12a contain a free pair [12] [43], so they belong to the type of eq. (B.6). According 
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to the rule above, we pull out the free factor and then we identify the harmonics pairwise. 
The identification 1 = 2, 3 = 4 annihilates the graphs, so they cannot contribute to A\. On 
the contrary, the graphs survive the identification 1 = 3, 2 = 4, so they contribute to A^- 
If we now draw the gluons, we see that these graphs are reduced to lie, llf rotated by 90°. 

Analyzing all the graphs along the same lines we eventually identify the coefficient 
functions they contribute to. The complete identification is indicated in Fig. 12. Drawing 
the gluon lines brings us to the diagrams in Fig. 13 for Aq and Fig. 14 for A\. 





Figure 13. Diagrams contributing to the Aq function for k = 4 case. 



A couple of comments are now in order. First of all we notice that the configurations of 
gluon lines appearing in Figs. 13 and 14 are the same as the ones we have already selected for 
k = 3. Moreover, a diagram with a given configuration of interacting lines can contribute 
to different functions according to the structure of free propagators which dress it. In fact, 
all the allowed configurations of interacting lines are present for both functions. 

Now, computing the contribution of each graph is simply a matter of combinatorics. At 
large N, neglecting the space-time structure, for the Aq function we obtain 

^(13) = £(13) = ^ N 8 

2C< 13 ) = 4 4 AT 8 

0,(13) = 44^8 

Cif = C4 13 ) = 2 x 4 4 iV 8 

2Cf 3) = 2 x 4 A N 8 (B.7) 
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whereas for A\ 

C (14) = c (14) = ^(14) = = 44^8 

2Ci 14 ) = 4 4 iV 8 
2C} 14) = 4 4 iV 8 

£,(14) = 44^8 

< ] = eg? = 2 x 4 4 iV 8 (B.8) 

It is important to notice that, at large N, the combinatorial factors only depend on the 
structure of the interacting lines: Diagrams with the same interactions but different config- 
urations of free propagators have the same coefficient. 




Figure 14. Diagrams contributing to the A\ function for k = 4 case. 
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